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Abstract: We investigate five-dimensional static (non-)extremal black hole solutions in 
higher derivative Anti-de Sitter gravity theories with neutral scalars non-minimally coupled to 
gauge fields. We explicitly identify the boundary counterterms to regularize the gravitational 
action and the stress tensor. We illustrate these results by applying the method of holographic 
renormalization to computing thermodynamical properties in several concrete examples. We 
also construct numerical extremal black hole solutions and discuss the attractor mechanism 
by using the entropy function formalism. 
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1. Introduction 

Among the fundamental interactions, gravity is very special. Gravity couples via a dimen- 
sional coupling constant, the Newton constant Gn, and so it is intrinsically non-renormalizable. 
Irrespective of the fundamental nature of quantum gravity, the gravitational low-energy de- 
grees of freedom are encoded in the metric of spacetime itself. However, there is no reason 
to believe that the effects of our present theory are the whole story at the highest energies. 
Indeed, non-renormalizability can be interpreted as a natural feature of a theory for which 
the action is not fundamental but arises as an effective action in some energy limit. At high 
enough energies — for sufficiently strong curvatures and sufficiently small distances — new 
interactions and new degrees of freedom will be required. 
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The fact that the gravitational action is proportional to R and only R is not due to any 
symmetry and, unlike other theories, can not be argued on the basis of renormalizability. 
Indeed, the low energy effective gravity action that obeys principle of equivalence and general 
covariance has a generic structure. That is the usual Einstein action plus a series of all possible 
interactions which are consistent with general covariance and local Lorentz invariance, i.e., 
higher curvature terms and also higher derivative terms involving the 'low-energy' matter 
fields. 

The effects of heavy particles appear to be local interactions when viewed at low energy. 
That is the fields at different spacetime points are independent degrees of freedom with inde- 
pendent quantum fluctuations. One important caveat related to the interpretation of gravity 
as a (local) effective field theory is as follows: in a local field theory one expects an entropy 
proportional to the volume, but that is not true for black holes. In classical gravity, a fixed 
energy-density in a sufficiently large volume will collapse into a black hole. 1 However, the 
holographic principle [1] was proposed to rescue this situation: gravity in D dimensions is 
equivalent with a local field theory in D — I dimensions. The AdS/CFT correspondence [2] 
(see [3] for a nice recent set of lectures) is a concrete realization of the holographic principle. 
Such correspondence is referred to as duality in the sense that the supergravity (closed string) 
description of D-branes and the field theory (open string) description are different formula- 
tions of the same physics. This way, the infrared (IR) divergences of quantum gravity in the 
bulk are equivalent to ultraviolet (UV) divergences of dual field theory living on the boundary. 
When we specify the CFT and say on which space it lives we are implicitly providing a set of 
counterterms for the gravity solution. These counterterms are local and depend only on the 
intrinsic boundary geometry [4, 5] (see, also, the reviews [6]) — a different method was pro- 
posed in [7]. Thus, one can compute the thermodynamical quantities in the gravitational side 
by employing the quasilocal formalism of Brown and York [8] supplemented by the boundary 
counterterms. The connection between the holographic charges and the various alternative 
definitions of conserved charges in AdS was explored in [9]. 

In studying string theories at low energy scales, the massive states may be integrated out 
to yield an effective action for the massless modes, with the same symmetries as the original 
string theory. Thus while the (super)gravity action is unique if we restrict to terms with two 
derivatives, interactions quadratic or higher order in the curvature tensor are allowed by the 
symmetries and so appear as well. However, such terms will require a dimensional constant 
to appear along with the derivatives. In string theory this constant turns out to be a', the 
inverse string tension. 2 

In this paper we investigate charged AdS black holes in the presence of higher derivative 
terms. We must note that, unlike in general relativity, in the presence of higher derivative 

1 In quantum gravity, the existence of local operators is problematic due to the causality. It is well known 
that the commutator of space-like separated local operators should be zero. However, since the gravity is 
dynamical the metric itself fluctuates and so the space-like intervals are not well defined. 

2 This constant defines what is meant by 'slowly varying fields' in the sense that the derivative corrections 
may be ignored for fields that are slowly varying on the scale of the string length l s ~ \/a'. 
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corrections there are two families of solutions. We propose counterterms that regularize 
the action and the stress tensor of both branches (for horizons with spherical, toroidal, and 
hyperbolic topologies). 

We obtain the stress tensor and the conserved charges for exact static charged non- 
extremal black hole solutions with Gauss-Bonet (GB) term and find perfect agreement with 
Wald formalism [10]. In the extremal limit we explicitly show that the near horizon geometry 
of the solution remains AdS2 x S 3 after including a'-corrections. The results we obtain provide 
a robust check of the entropy function formalism [11]. Indeed, we find that, for our exact 
solutions, the radius of AdS2 receives corrections but the near horizon geometry remains 
AdS2 x S 3 (H 3 ). In this way we obtain the generalization of Bertotti- Robinson geometries 
[12] with GB term. 

We also apply the counterterm method to 5-dimensional charged black hole solutions in 
gravity theories with £7(1) gauge fields and neutral scalars. We obtain numerical solutions 
and generalize the results of [13] by including the higher derivative terms. 

In the extremal limit we study the attractor mechanism by using the entropy function 
formalism [11, 14]. This method is based on the near horizon geometry and its enhanced 
symmetries but does not provide a proof for the existence of a complete solution in the bulk. 
For some special values of the couplings, we present numerical solutions with a finite horizon 
- this confirms the results in [15] where the equations of motion in the bulk were solved 
perturbatively order by order. Thus, we can safely apply the entropy function formalism. 

An overview of the paper is as follows: in section 2 we study in detail the AdS charged 
black holes with GB term. We compute the stress tensor and the conserved charges of the 
exact non-extremal black hole solution by using the counterterm method and compare with 
the results obtained by Wald formalism. We present a preliminary discussion on the thermo- 
dynamics in both, canonical and grand-canonical ensembles. We also study the extremal limit 
and interpret our results within the entropy function formalism. Section 3 is dedicated to 
studying black hole solutions in AdS gravity with U(l) gauge fields non-minimally coupled to 
scalars in the presence of GB term. We present numerical non-extremal solutions and discuss 
in detail their properties by using the counterterms proposed in section 2. In section 4 we 
study the extremal limit in the case of massless scalar fields, construct numerical solutions, 
and investigate the attractor mechanism for these solutions. In section 5, we discuss our 
results. An appendix gives some calculational details on Wald formalism. 

2. Charged AdS black holes with Gauss-Bonet term 

In this section we compute the conserved charges of AdS charged black holes with GB term 
by using both, the counterterm method and Wald formalism. The GB term is a very natural 
correction term to the Einstein action in the sense that the equations of motion contain no 
more than second derivatives in time. The main reason we are interested in GB term is due 
to the existence of exact solutions. In the extremal limit, we explictly check that the near 
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horizon geometry still remains AdS2 x S 3 (H 3 ). We also use the entropy function formalism 
to interpret our results. 



2.1 Non-extremal case 



AdS spacetime is a maximally symmetric Lorentzian space (i.e. the number of Killing vector 
fields is the same as for flat spacetime) with constant negative curvature — in D-dimensions 
the symmetry group is SO(D — 1, 2) and the topology is AdSo = R D ~ 1 x S l . AdS spacetime 
arises as the natural ground state of gauged supergravity theories and plays an important 
role in understanding holography in string theory. 

As we are interested in AdS gravity with higher derivatives (see [16] for a recent review), 
we begin by establishing our conventions for the action. Exact solutions are presented in [17]. 
Note, however, that our conventions differ from the ones in [17] and we also correct some 
important typos. 

In this section, we will focus on a 5-dimensional theory of gravity with negative cosmo- 
logical constant coupled to vector field, whose general action has the form 



/ = - -L I d 5 x^[R -2A- F^ u + a'L GB ] 

K 5 JM 



(2.1) 



where = d^A u — d v A^ is the gauge field, A is the cosmological constant, and K 2 = \6irG. 
We use Gaussian units so that factors of 4tt in the gauge fields can be avoided. The GB term 
L GB = R 2 - AR^R^ + R^BF"*? appears in the low-energy effective bosonic string theory 
— in type IIB superstring the leading corrections are cubic in a'. 

Within this theory there is a straightforward generalization of the Reissner-Nordstrom 
(RN) solution 



with 



ds 2 = -N(r)dt 2 + Nir^dr 2 + r 2 dY>\ 



(2.2) 
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(2.3) 



where $ is a constant which is chosen such that A t (rh) = and is the largest positive 
root of N(r) that is typically associated to the outer horizon of a black hole — note that 
the condition N'{rh) > implies the existence of a minimal allowed value of r^. Here, L is 
the radius of AdS spacetime and it is related to the cosmological constant by A = — 6/L 2 
and k = 1,0,-1 corresponds to black holes with spherical, planar, and hyperbolic horizon 
topologies. The expression of N(r) has an extra parameter e = ±1, that implies the existence 
of two branches of solutions. 

Let us discuss now some known limits of the solutions (2.3) — more details can be 
find in [17]. The minus-branch solution reduces in the limit of a' — > to the RN solution 
of the Einstein-Maxwell-A system, i.e. N(r) = k — m/r 2 + q 2 /r 4 . On the other hand, 
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diverges for the e = +1 branch and so there is no smooth limit in this case, since N(r) = 
r 2 /(2a') + k — m/r 2 + q 2 /r 4 - r 2 /L 2 as a' -> 0. 

The background approached asymptotically by these solutions corresponds to an AdS§ 
spacetime with an effective radius 



L e// = Lyi±^, where U=Jl~ (2.4) 

This limit (m = q = 0) corresponds to AdS with higher derivatives. This effective radius 
of AdS with higher derivative corrections will play an important role in the subsection 2.1.1 
where we will define the counterterms for the action and the stress-energy regularization. It 
is clear that by adding higher derivative corrections (even for small a') the theory contains 
new solutions (in our case a new branch) unavailable in general relativity. 

2.1.1 The counterterm method 

We start by reviewing some known useful facts about the quasilocal formalism of Brown 
and York [8]. The gravitational field (the metric tensor) couples to the energy momentum- 
tensor (or stress tensor) of every other field in nature. In general relativity, mass is merely 
one aspect of the stress tensor, and gravitational energy is non-local as follows from the 
equivalence principle. By choosing a coordinate system that is inertial in a given volume 
element one can make the stress-tensor vanishing (since T l kl is vanishing). Thus, it has no 
meaning to speak of a definite localization of the energy of the gravitational field in space. 
One can measure the gravitational field by the geodesic deviation of two observers — a single 
observer cannot distinguish it from kinematical effects. In other words, curvature cannot be 
measured on a point line, but requires a 2-surface at least. Since an appropriate definition of 
the gravitational energy cannot be found locally, a quasilocal definition is sought. 

One way to compute the energy of a gravitational system is by enclosing it with a surface 
- the observers living on this surface can make measurements and compare the results. In 
the quasilocal formalism, the surface stress tensor for spacetime and matter is defined by 

where 7 a b is the induced metric on the enclosing surface. Even if there are similarities be- 
tween the definiton of the matter stress tensor and the boundary stress tensor, it is worth 
to emphasize that T ab characterizes the entire system, including contributions from both the 
gravitational field and the matter fields [8] . 

As usual in gravity theories, the action (2.1) should be supplemented with suitable bound- 
ary terms to obtain a well-defined variational principle. For Einstein gravity, one considers 
the Gibbons-Hawking surface term [18] 



r(E) 

* 1 



b 8vrG 



l — [ d 4 x^K (2.6) 
TLr JdM 



- 5 - 



where 7^ and K are the induced metric and the trace of the extrinsic curvature of the 
boundary, respectively. 

A similar term occurs for Gauss-Bonnet gravity and reads [19] 

where is the four-dimensional Einstein tensor of the metric ^ ab and J is the trace of 

J ab = i (2KK ac K c b + K cd K cd K ab - 2K ac K cd K db - K 2 K ab ) (2.8) 

Variation of the action I + I b + I b now gives an expression that does not contain normal 
derivatives of 5g ab . 

It is well known that the total action contains divergences even at tree-level — they 
arise from integrating over the infinite volume of spacetime. We regularize the divergences by 
using the procedure proposed in [5]. This technique was inspired by the AdS/CFT duality 
and consists in adding suitable counterterms Id to the action of the theory in order to ensure 
its finiteness. 

We have found that the action of the solutions in this paper can be regularized by the 
following counterterm 

I « = ^[ d 4 xv^( Cl -|R) (2.9) 
o^G J dM 2 

where R is the curvature scalar associated with the induced metric 7. The consistency of the 
procedure requires 

Cl = --^-(2 + eU) C2 = ^£l(2-eU) (2.10) 

L eff 2 

This counterterm is general and can be used to regularize the action of both branches. 3 
For solutions with a well defined Einstein gravity limit, one finds as a' ^ 0, that c\ — ► 
-3/L + a'/L 3 + O(a') 2 , c 2 -► L/2 + 3a'/2L + 0(a') 2 that match the results in [21, 22] (see, 
also, [23]). 

Gravitational thermodynamics is then formulated via the Euclidean path integral, where 
one integrates over all metrics and matter fields between some given initial and final Euclidean 
hypersurfaces. Semiclassically the total action is evaluated from the classical solution to the 
field equations. The thermodynamical system has a constant temperature 

where (3 is the periodicity of the Euclidean time determined by requiring the Euclidean section 
be free of conical singularities. 



3 A discussion of the counterterm method for GB gravity with cosmological constant also appears in the 
forth-coming paper [20]. 



- 6 - 



To evaluate the action, one express the bulk action as a total derivative 

±(R - 2A + a'L GB - F 2 ) = ^~r 3 N' -^+6r a' (N — k) N'^j (2.12) 

where a prime for a metric function denotes a derivative with respect to the radial coordinate 
r. After adding the boundary terms, one finds that the Euclidean action is finite and contains 
two terms, I = I as + I eh . These two terms represent the contributions from the boundary 
and the event horizon and their expressions are 

I as = + C with I°> = k 2 3 ^ f J k (Sell - 2) (2.13) 

with Vk the area of the surface One can easily verify that the action computed according 
to a background subtraction coincides with the above expression up to the Casimir term Iq s 
(the background choice in this case corresponds to a q = vacuum EGB-AdS solution). 

Varying the total action (that contains the boundary terms (2. 6), (2. 7), and (2.9)) with 
respect to the boundary metric h a b, we compute the divergence- free boundary stress-tensor 



Tab — 



8ttG 
where 



(^Kab ~ K~j ab + dj ab + C 2 G a6 + ^(Qab ~ ^Qlabi^j (2.15) 



Q ab = 1KK ac K c h - 2K ac K cd K db + K ab {K cd K cd - K 2 ) (2.16) 
+2KR ab + RK ab - 2K cd R cadb - 4R ac K c b 

with R a bcd an d R a b denoting the Riemann and Ricci tensors of the boundary metric. 

Provided the boundary geometry has an isometry generated by a Killing vector £ l , a 
conserved charge 

5 = j> d 3 S i gTij (2.17) 

can be associated with a closed surface £ [5]. Physically, this means that a collection of 
observers on the hypersurface whose metric is hij all observe the same value of £2^ provided 
this surface has an isometry generated by £. The mass/energy M is the conserved charge 
associated with the Killing vector £ = d/dt. For charged black holes, the expression of the 
nonvanishing components of the boundary stress tensor are 

SttGT™ = (^mL eff - k 2 ^{2 - MJ^j ± + 0(l/r 6 ) (2.18) 
SnGTt = [~mL eff + 3k 2 ^(2 - 3eU)) 1 + 0(l/r 6 ) (2.19) 
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where w denotes an angular direction on S3 (note that in 1/r 4 order, this is a traceless stress 
tensor). 

Due to its high degree of symmetry, AdS space has a simple form in a large number 
of coordinate systems. By choosing different foliations of the spacetime one can describe 
boundaries that have different topologies and geometries (metrics), affording study of the 
CFT on different backgrounds. Specifically, we found additional Casimir-type contributions 
to the total energy depending on the slicing topology in accord with the expectations from 
quantum field theory in curved space. This can be seen for the solutions discussed in this 
section, whose mass computed according to (2.17) is 

OTA 

M = T^b' n + k2 e^ lV3cU - 2) < 2 ' 20) 

where the last term is the Casimir energy. 

The metric on which the boundary CFT is defined is found by getting rid of the divergent 
conformal factor, h a b = limj—xxj -^r-^f a b, and corresponds to 

h ab dx a dx b = -dt 2 + L 2 eff dT,l (2.21) 

If such a CFT exists the theory lies in the landscape of string theory and the bulk theory is 
manifestly consistent as an effective theory, otherwise the theory is part of the swampland 
[24]. 

2.1.2 Wald formalism 

One way of understanding black hole entropy comes from the use of Euclidean analog of a 
black hole spacetime. Whenever it is not possible to foliate the Euclidean section of a given 
(stationary) spacetime by a family of surfaces of constant time, gravitational entropy will 
emerge. Another approach to gravitational entropy is the Noether charge formalism of Wald. 
The relation between the two methods was explored in [25] (see, also, [26, 27]). 

When we add R 2 corrections to the action the entropy is no longer given by the area 
law — instead, to computing the entropy of the black holes (2.3), we will use a more general 
formula proposed by Wald [10]: 



= ~ 27T / d3x 7T5 e ^cd (2-22) 



m 

where 7i is the bifurcate horizon and e^ u is the binormal to the bifurcation surface. Inter- 
estingly enough, the entropy can still be expressed as a local functional evaluated at the 
(bifurcate) horizon. In this construction, the entropy was obtained from the Noether charge 
that is the integral of a 3-form associated with the diffeomorphism invariance of the theory. 
It is worth noticing that Wald formalism can be applied to non-extremal black hole solutions 
in generally covariant theories of gravity. 

The most general formula for the entropy for a Lagrangean of the form 

R 



-I 



d°x^/^g 



16vrG 



2A + aR l + PR^R^ + iR^ P aR^ pa 



(2.23) 
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is given by (see apendix) 

S = f d 3 xVh \l + 2K 5 aR + K b f3{R - h ij R tj ) + 2K^{R - 2h ij R ij + h^h kl R ikjl ) 



We are interested in GB term for which the expression for entropy becomes, 

d 3 x^h \l + 2a'h i ' k h j ' l R ijH 



in 
where 



(2.24) 
(2.25) 



a 



4a' 



for GB term and K$ = 16itG. It is easy to show that 



h ik hi l R ijkl = 4- (2.27) 
r + 

and hence entropy becomes 

S=^r h (rl + 12ka') (2.28) 
that matches with (2.32) obtained by the counterterm method. 

2.2 The grand canonical and canonical ensembles 

The results above make possible a discussion of the thermodynamic properties of these charged 
black hole solutions. In a very basic sense, gravitational entropy can be regarded as arising 
from the Gibbs-Duhem relation applied to the path-integral formulation of quantum gravity, 
which in the semiclassical limit yields a relationship between gravitational entropy and other 
relevant thermodynamic quantities. In this approach, the expression of the entropy is 

S = P{M- - I (2.29) 

upon application of the Gibbs-Duhem relation to the partition function, with chemical po- 
tentials Cj and conserved charges /Xj. For the situation in this work, £ corresponds to the 
electrostatic potential <!>, while \i is the electric charge Q, with 

$ = ^T' ( 3 = i%2v / 3g (2.30) 

To compute the entropy, it is convenient to express everything in terms of (rh,q) 

TH ~ 2vrL2 r 2 /L 2 + W/L 2 > m- ^+kr h + L2 +2ak (2.31) 

the action being given by the sum of (2.13) and (2.14). In this way, one finds the following 
expression for the black hole entropy: 

S=^r h (r 2 h + 12ka>) (2.32) 
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One can easily verify that the first law of thermodynamics dM = TndS + QdQ also holds. 

The corresponding equation of state (analogous to f(p,V,T), for, say, a gas at pressure 
p and volume V) reads 

rp 1 /Q 3Q + 2c&(3fc-4<l> 2 )L 2 

H 6vrL 2 V$ Q+16fca'$ " 1 j 

A discussion of the corresponding thermodynamical properties can also be approached. In a 
grand canonical ensemble one finds the Gibbs free energy 

W[T H , $] = M - T H S - Q$ = W + Wi, where W = -^r — f L(L + V^ 2 - 8a') - 8a') , 

SttG 16 V / 

Vfc Q 2 (3Q + 4$L 2 (4$ 2 - 3fc)) + 48a'fc$Q(9Q + 4$(3fc - 4$ 2 )L 2 ) 
1_ 8vrG 96L 2 $ 2 (Q + 16fca'$) ^ j 

where Q is given as Q(Th, 3>) by the equation of state (2.33). 

One can consider insted a canonical ensemble, where the temperature and electric charge 
are keept fixed. The Helmholtz potential F = M — TS in this case is 

Vfc 3A: 2 



F[T H ,Q]=F + F 1 , F = ^— [L(L + Vl 2 - 8a') - 8a') , (2.35) 

Vfc Q 2 (-3Q + 4$(3fc + 20$ 2 )L 2 - 144a'/e$Q(3Q + 4$(/e - 4$ 2 )L 2 ) 
1 ~ 8^G 96L 2 $ 2 (Q + 16fca'$) 

where the electrostatic potential $ is given as a function of Th,Q by the equation of state 
(2.33). 



2.3 Extremal case 
2.3.1 Exact solution 

The non-extremal black holes have a non-zero temperature that can be evaluated by eliminat- 
ing the conical singularity in the Euclidean section. Once we impose the periodicity condition, 
the Euclidean time circle closes off smoothly and the Euclidean geometry becomes a 'cigar'. 
On the other hand, an extremal Euclidean black hole has a different topology. That is an 
infinite long throat for which the Euclidean time circle does not close off. In this case, one 
is forced to work with an arbitrary periodicity of the Euclidean time leading to ambiguos 
results (though, see [28]). 

However, on the Lorentzian section the picture is quite satisfactory: an extremal black 
hole is obtained by continuosly sending the temperature of a non-extremal black hole to zero. 
While the temperature vanishes, the area of the horizon can remain finite. 

Thus, to obtain the extremal black hole solution we work on the Lorentzian section. The 
extremal limit can be equivalently obtained by imposing the constraint that the horizon is 
degenerate (i.e., N(r) has a double root: N(rn) = N'{th) = 0). 
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One can easily solve the equations system to obtain: 

m = 2 (*(te' + r |) + |^ (2.36) 

«•-"*(£ + £) (2.37) 

where rjj is the horizon radius. 

Using the method of [29] (see section 4.2) it is straightforward to show that the near 
horizon geometry of (2.3) is AdS-z x S 3 (H 3 ) and just the radius of AdSz receives corrections 
(we will compute and provide its concrete value in section 2.3.2). 

These geometries are interesting in their own right and provide the generalizations of 
Bertotti-Robinson geometries with GB term. These solutions are the topological product of 
two manifolds of constant curvature. They are conformally flat and are supported by a flux 
through S 3 {H 3 ): 



ds * = -F(r)dt> + ^ + RldXl with Hr) = k + 4 ^ + \ k k ^ (2.38) 

and satisfy the equations of motion with a gauge field At = r^J\2/L 2 + 6kjR\ — here is 
a constant and so the size of E| is fixed. Since H 3 is not compact, the solution also exists for 
a vanishig electric potential in this case. 

2.3.2 Entropy function 

Wald's construction for the entropy can be used in any general coordinate invariant theory 
of gravity including those with higher derivative terms in the action. The black hole entropy 
is obtained in terms of the field configurations near the horizon — all the physical fields (not 
just the curvature) and their derivatives should be regular at the bifurcation surface. This 
method is based on a Lagrangean derivation of the first law of black hole thermodynamics and 
can be used directly on the Lorentzian section (no 'Euclideanization' is required). However, 
this method can be applied just to black holes with bifurcate Killing horizons and so not to 
extremal black holes. 4 

This method was also extended by Sen to extremal black holes and it is reffered to 
as the entropy function formalism [ll]. 5 Extremising the entropy function is equivalent to 
the equations of motion in the near horizon limit and its extremal value corresponds to the 

4 It is also worth emphasizing that a black hole formed by a collapse process does not have a bifurcate 
horizon. 

5 It is known that the near horizon geometry of stationary extremal black holes contains an AdS2 space 
[30] (see, also, [31]). The entropy function is constructed, on an SO(2, 1) x SO(3) (for static black holes) or 
SO(2, 1) x U(l) (for stationary black holes) symmetric (near horizon) background, by taking the Legendre 
transform (with respect to the electric charges and angular momentum) of the reduced Lagrangian evaluated 
at the horizon. 
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entropy. A discussion on the entropy function formalism and the Euclidean section method 
can be found in [32]. 

In this section we use the entropy function formalism and compare the results with the 
ones in the previous subsection. The general metric of AdS 2 x S 3 can be written as 



ds 2 = v 1 (-p 2 dT 2 + — dp 2 ) + v 2 dVL\ , 



(2.39) 



The field strength ansatz is F = edr A dp and, for this geometry, the GB term comes out to 
be GB= —2A/v\v 2 . Thus, the entropy function F(vi, v 2 , e, Q) is given by 



F(vi,v 2 , e, Q) = 2ir[Qe - f(v 1 ,v 2 , e)] , 



(2.40) 



f(v 1 ,v 2 ,e) = 2tt 2 
The attractor equations are: 



2v^ 2 + Gvit/v^ + 2— — e l + v\v\ 



3/2 



3/2 



(j^j - 24a' ^T 2 



9F . 2 2 , 2 /12 



<; „ 6v{ - 2v 2 e 2 + v{v 2 1-^1 =0, (2.41) 

- n _ ...... , ..2 , „._„2 , U?«2 f 12' 

<9t> 2 



= ^ -vi«2 + v\ + v 2 e 2 + ^ -4a' Vl = 0, (2.42) 



flT7 3/2 

££ = Q = 87r 2 ^e. (2.43) 
oe v\ 

Let us now discuss in detail these equations. One important observation is that by 
adding the GB term to the action just the second attractor equation is modified. One can 
easily eliminate v\ from the first equation by using the third one and so the value of v 2 does 
not change — we obtain the following relation between the electric charge and the horizon 
radius (y 2 = r 2 H ): 

Using the conventions from the section (2.1) and the relation between the physical electric 
charge and the charge parameter Q = \^3q we can see that this relation matches (2.36). It is 
worth emphasizing that just the radius of AdS 2 receives corrections: 

Vl = 4TT3V^ ( } 

Replacing v±(v 2 ) and e(Q,v 2 ) back in the entropy function we obtain the entropy of the 
extremal black hole: 

, , / 12a' \ 

Sextremal = S^rjj ^1 + -5- J (2.46) 
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The entropy of extremal black hole has the same form as for the non-extremal black hole 
(2.32) (Vfe = 2-7T 2 and = l/167r), though the radius of the horizon (r#) is different. Not 
surprisingly, the form of the a' correction is the same — both methods, Noether charge and 
entropy function, are based on a Lagrangean derivation. 

In general there are two types of first-order corrections due to higher derivative terms. 
The entropy/area law is modified due to the additional terms in the action 6 and/or the 
modification of the area due to the change of the metric on the horizon (the extra terms 
in the action may change the equations of motion). In our case, the horizon radius, V2, of 
the extremal black hole remains unchanged after adding the GB term and so the entropy is 
changed due to the suplementary terms in the action. 

3. Charged black holes with scalar hair 

In this section we generalize the results of [13] by including the GB term. We obtain numerical 
solutions 7 and discuss how the GB term affects their properties by using the counterterms 
proposed in the previous section. 

3.1 The model 

We consider the generalization of the RN black holes in a five-dimensional theory of gravity 
coupled to a set of scalars and vector fields, whose general action has the form 

I[G^, Ajf] = J d 5 x^[R + a'L GB - Gu^d^ 1 d»<j> J 

-U^F^^-Vm, (3.1) 

where F£ v with A = (0, • • • N) are the gauge fields, 4> = (4> L ) with (I = 1, • • • , n) are the 
scalar fields, V(4> t ) is the scalar fields potential. 

The equations of motion for the metric, scalars, and the gauge fields are given by [17, 15] 

iV - X -R 9ilv + g^V^) + a'H^ = 8vrGT™ a " er (3.2) 

(3.3) 

d,[V^(fAB(4>)F B ^)} =0 (3.4) 
where TJ^ atter is the matter stress tensor and H^ v is given by [17] 

= 2(RR flu - 2R^ a R" - 2R a ^R mv p + R^f 1 R va i3-i) - -^S^uLgb (3.5) 

6 These terms are evaluated using the zeroth order solutions for the metric and the other fields 
7 We thank Eugen Radu for advice in finding the numerical solutions and explaining to us the methods used 
in [33]. 
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The Bianchi identities for the gauge fields are F A r.^ = 0. 

We assume that the scalar fields approache asymptotically constant values, (j) 1 ^, which 
corresponds to an extremum of the potential such that dV/d(j)\^ = and V((/>oo) = — 12/ L 2 < 
0, with the expansion 



1 d 2 V 

v(4>) = v( ( p OD ) + -- 



2 d^d(f> J 



+ ... (3.6) 



the scalar field masses being set by d 2 V '/d^ 1 d(f> J \^ = hij. 

Under these asumptions, the background of the theory is given by the solution 

ds 2 = -(k + -^)dt 2 + — + r 2 dS| (3.7) 

L eff k+j^r- 

L eff 



with 4> = 0oo, k = ±1,0 and the effective length scale is L e ff = L^J\ + U /2, with U = 
y/l - 8a' /L 2 as in the case of EGB-A theory. 

Restricting to static solutions, we consider the metric ansatz 

dr 2 

ds2 = W^) + ^ " N ^< r f dt ^ ( 3 - 8 ) 

and a purely-electric abelian field ansatz A B = W B (r)dt, the scalar fields being also functions 
only of the radial coordinate r. 

The Maxwell equation implies the existence of the first integrals 

W' A = f AB q B ^ (3.9) 

where qs are constants fixing the electric charges of solutions, Qb = qBVk/(2irG) and f AB is 
the inverse of Jab- The electric potentials are the integrals of F B , being fixed up to arbitrary 
constants Q B which are chosen such that Af vanish on the event horizon. 

It is more convenient to combine the equations of motion (see [15]) to obtain the following 
equivalent system of differential equations: 

l r (rN' + 2N- 2k) - 3a'(N - k)N' + ^NGutf 1 '(f>' J + ^r 3 F(0) + ^f AB q A q B = 
a' = 1 rHhjf^_ 

a 12 r 2 /4 - a'(N - k) 1 j 

The first equation does not contain any second derivatives and is the Hamiltonian con- 
straint. We notice that the equations of motion can also be derived from the one-dimensional 
reduced Lagrangean: 



L red = -3ra(rN' + 2N - 2k) + 12a' (iV - k)aN' - r^aNGutf V - r 3 aV((f>) - ^f AB q A qB 
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We are interested in black hole solutions approaching asymptotically the background (3.7), 
that suggests to use the following form of the metric function, N(r): 



N(r) = k-^ + ^- (3.11) 
The first equation in (3.10) implies that m(r) satisfies the following equation 
\ (Um + Ll ff (U - = ViVG^V' + \AV(<P) - V(M) + 

The computation of the action and the stress tensor of these configurations can be done 
by using a similar approach to the one discussed in section 2 and we will not present the 
details here (without a' corrections, see [13]). Similar to the case without scalars, the volume 
term in the action (3.1) has a total derivative structure and so it can be expressed in terms 
of the difference of two surface integrals. 

The divergencies associated with the asymptotic AdS structure of the solutions can be 
removed by supplementing (3.1) with the same boundary counterterms (2.9) as in section 2. 
After Wick rotating t — > it to the Euclidean section, we found that the action can be written 
in the usual 'quantum statistical' form 

I = P(M- Q A $ A ) ~ J^r h {r 2 h + 12a'k) (3.12) 

where (3 is the periodicity of the Euclidean time. Here, M is the mass of the black hole and 
we will present its values for some concrete examples in the next sections. 

The value of (3 is arbitrary for soliton solutions (that exist in the absence of gauge fields, 
e.g. [34]) or extremal black holes. However, the regularity of the Euclideanized solutions as 
r — > rh imposes 

P = h = N>{r%{r h ) (3 - 13) 

for non-extremal black hole solutions. 

The mass of these solutions can be computed within the quasilocal formalism by using 
the generic relation (2.17), where the boundary stress tensor is still given by (2.15). However, 
the situation is different for theories with massless scalar fields and in the presence of massive 
scalar fields — we shall discuss these cases separately. 

The explicit construction of solutions requires specification of the functions Gjj and /ab- 
In what follows we consider a model with one single scalar (i.e. Gjj = 25udu, 4> l = <fi) and 
two gauge fields with modulus dependent couplings of the form 

f AB {4>) = S AB e aB ^ (3.14) 
Moreover, we shall restrict to solutions with a smooth Einstein gravity limit. 
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3.2 Unattractor solutions with a massless scalar field 

In their simplest version, these solutions have a constant value of the scalar potential, V(4>) = 
2 A = — 12/L 2 , which is the case considered here. The generic solutions have a non-degenerate 
horizon and are easier to study. Near the event horizon, they admit a power series expansion 
of the form (here we restrict to the first terms in the series) 

N(r) = f 1 (r-r h ) + ..., a(r) = a h + ^f^* (r - r h ) + . . . , (3.15) 

<H0 = <t>h - Tri^ai^N? + a 2 e- Q2< N|)(r - r h ) 
2r hh 

where 

7l " 3^ + 4^) (3 - 16) 

The coefficients of all higher order terms in the expression of N, a, (ft are fixed by the two 
parameters 4>h,&h- One can easily see that the condition f\ > imposes the existence of a 
minimal value of for given values of <ph,Qi,Q2- 

One can also construct an approximate solution at the boundary in terms of three free 
parameters (poo, S, and Mq 

M(j 2 e - ai 0oo 2 + e -a 2 ^ 2 

iY(r) = 1 -^ + zj; + sZT^ + •••' (3 - 17) 

<A(r) = <Aoo + ^ + ..., ^) = 1 "^8+--- ( 3 - 18 ) 

The next leading term at the boundary corresponds to a normalizable mode and the 
black hole is a state in the boundary CFT. 

By applying the quasilocal formalism discussed above, one finds the mass of these solu- 
tions 

while the entropy of solutions has the same form as in the RN case, S = + 12a' k) 

We will explictly check by our numerical analysis that M, rh, and 4>h depend of the 

asymptotic boundary data (<poo)- This is in contrast with the extremal case (see section 4) 

where we obtain an attractor behaviour of the horizon. 

Although an exact solution of the equations of motion (3.10) appears to be intractable, 

here we present arguments for the existence of non-trivial solutions which smoothly interpolate 

between the asymptotic expansions (3.15) and (3.17). 

Starting from the event horizon expansion (3.15) we integrated the equations towards 

r — > oo. The integration stops when the asymptotic limit (3.17) is reached with a reasonable 
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accuracy. In this approach, the input parameters are k,rh,qi,q2,oii,a2,L and the value 
(f>h of the scalar field on the horizon. The equation for a decouples from the rest, and the 
requirement that a — > 1 as r — > oo can be relaxed during the numerical integration, cr(rh) 
subsequently being multiplied by an appropriate constant factor so that the correct asymptotic 
behaviour is recovered. 

We have solved the equations of motion for several values of a,\ = — «2 = 2a and a 
large set of rh,qi,q2,4'h,L. We follow the usual approach and, by using a standard ordinary 
differential equation solver, we evaluate the initial conditions (3.15) at r = + 10 -4 for 
global tolerance 10~ 12 , for a fixed parameter (p^ and integrating towards r — > oo. 

The complete classification of the solutions in the space of parameters is a considerable 
task that is not aimed in this paper. Also, we shall restrict to the case of spherical topology 
horizon, although we have found topological black holes as well. 
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Figure 1: The profiles of the functions m(r), a(r) and cf>(r) are shown for typical k = 1 non-extremal 
black holes with with rt = 1, qi = 0.3, q2 = 0.5, A = —6, a' = 0.1 and two values of the scalar field on 
the event horizon 

For all configurations we have studied, the metric functions N(r), cr(r), and the scalar 
(j){r) interpolate monotonically between the corresponding values at r = and the asymptotic 
values at infinity, without presenting any local extrema. In Figure 1 we plot the profiles of 
two typical solutions with different values of 4>{rh)- The evolution of the solution data as 
a function of the event horizon radius is reported on Figure 2. For small values of the 
numerical analysis strongly suggests that an extremal black hole solution is approached for a 
critical non-zero value of the event horizon radius. 

3.3 Configurations with a massive scalar field 

The situation is different when allowing the scalar field to present a mass term. As discussed 
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Figure 2: The relevant parameters are ploted as a function of the even horizon radius for k = 1 
non-extremal black holes with A = —0.001, a' = 0.01, qi = 0.3, q2 = 0.5, a\ = — a-2 = 2. 

in the last years by various authors, the field equations (3.10) with a' = admits a variety 
of solutions. This includes also configurations where the asymptotic behaviour of the scalar 
field is assumed to be slower than that of a localized distribution of matter. By relaxing the 
standard asymptotic conditions for asymptotically AdS solutions, it is possible to preserve 
the original symmetries at infinity, while the conserved charges are modified by including 
matter field terms (see e.g. [35] and the references there). 

It would be interesting to see how these features are affected by the presence of a GB 
term in the action. In this context, we start by discussing the issue of Breitenlohner-Freedman 
(BF) bound for a scalar field in the presence of a GB term in the Lagrangean. One can easily 
see that for a' > 0, the mass bound increases according to [i 2 BF = —A/L 2 ^. For example, 
the generic asymptotic behaviour of the scalar field (j){r) in the background (3.7) is 

where 0i, 4>2 are constants and 

A± = 2 M ± ^/i-mV/^bf) ( 3 - 21 ) 

Here we assume ^i 2 L^^ + A > 0. Imposing that both the A_ and A+ solutions be normahzable 
results in a supplementary conditions on the parameter ix 2 L 2 ^+2> < 0. For fields that saturate 
the BF bound, A+ = A_ and the solution is <p(r) = 4>i/r x + log r/r A . 

We shall restrict here to a scalar field (j) that is tachyonic (fi 2 < 0) and its mass is in the 
BF range, i.e. decaying at infinity according to (3.20) and not saturating the BF bound. We 
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also consider solutions with a well defined Einstein gravity limit. 8 

In this case, the conserved charges are well defined and finite despite the fact that the 
scalar field falls off slower than usual. The analysis in the procedure of holographic renormal- 
ization is based on finding the most general asymptotic solution of the field equations [37] and 
so the fall off of the matter fields is not an input in the computation. 9 Considering now the 
question of action and total mass-energy of these solutions, one can see that for 02/0, due 
to the back reaction of the scalar field, the boundary counterterms (2.9) are not enough to 
cancel all divergences in the on-shell action of the solutions. Therefore one has to supplement 
I c t with matter counterterms [38]. The counterterms (which include terms beyond the purely 
gravitational ones) are local and the action (and, also, the holographic stress energy tensor) 
has a contribution from the scalar fields. 

As a result, the total action of the solutions is 

-±(rl + 12ka'r h )V k (3.22) 

where Mq is the usual mass parameter, = A-^i^/^e// ( see ' e -§-' [36]), and the Casimir 
contribution is 

M casimir = k 2 ^^V k (3U - 2) (3.23) 
b47rG 

Thus the mass contains a supplementary term due to the slower decay of the scalar field. 
This dependence is crucial for the stress tensor to satisfy the correct Ward identities and to 
lead to sensible results. 

However, from the Gibbs-Duhem relation (2.29) one finds that the entropy of these solu- 
tions is still given by the relation (2.32). 



1 = ^(^ + M )+M mir -Q 1 <D 1 -Q 2 <D 2 



4. Extremal solutions and attractor mechanism 

The numerical extremal solutions can also be discussed by using similar methods as in previous 
section. However, the situation in this case is more involved. Similar to the non-extremal 
case, one may write an approximate form of these configurations near the event horizon — for 
static extremal black holes the near horizon geometry is AdS^ x S 3 . Let us first investigate 
the near horizon geometry of these black holes by using the entropy function formalism. For 
simplicity, we are again considering a theory with one scalar field and two U(l) (electric) 
gauge fields with the couplings given by (3.14). The general metric of AdS 2 x S 3 can be 
written as 

ds 2 = Vl {-p 2 dT 2 + ^dp 2 ) + v 2 dtil (4.1) 

8 A detailed discussion on the asymptotic form of the metric and the implications for the no4iair theorem 
can be find in [36]. In our case, the only difference is that the solutions approach at the boundary an AdS 
space with a modified radius, L e ff given by (2.4). 

9 We thank Kostas Skenderis for discussions on holographic renormalization method for this case. 
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The field strength ansatz is F A = e A drf\dp. Thus, the entropy function F{y\,v 2 , e A , qA, 4>h) 
is similar with the one in section 2, except that we have now non-trivial couplings between 
scalars and the U(l) fields: 



F(vi,v 2 ,e A ,q A ,4>h 
f(v 1 ,v 2 ,e) = 2vr 2 



= 2n[q A e A - f(v u V2, e A , q A , </> h )] , (4.2) 

3/2 / 1 9 \ 

-2vl 12 + Q Vl ^T 2 + <l V -2-f AB {4> h )e A e B + Vl v 3 2 /2 (72 J - 24a' ^ 



By using the same trick as in section 2 we can compute the horizon radius and the value 
of the scalar at the horizon by solving the following equations 

r 2 

a ie - ai<ph qi + a 2 e _Q2 NI = 0, e~ a ^ h q\ + e'^qj = 3rf,{k + -|) 

in terms of qi,q 2 and a±, a 2 . 

While 4>h = a2 ^L a \og(—((*2q 2 ) / (&iq 2 )) , the expression of Th(qi, q 2 , L) is very complicated 
and we do not present it here — note though that the relation between the horizon radius and 
the charges is similar with (2.44) where Q 2 is replaced by q\q 2 . The horizon value of the scalar 
does not depend of the boundary value (j)^ and so the near horizon geometry is universal. 
Consequently, the entropy of the extremal black hole does not depend of the boundary values 
of the scalar field. The behaviour of the scalar field is illustrated in figure 3 and the attractor 
mechanism is a direct consequence of the extremality condition. 

At this point, it is worth trying to find a whole solution interpolating between the horizon 
and the boundary — the entropy function assumes the existence of such a solution but does 
not prove it. 

Unlike in the non-extremal case, the value 4>h of the scalar field on the horizon and the 
event horizon radius are not free parameters and so the horizon data contain two essential 
parameters. The leading terms in this expansion read 



N(r) 



A(kL 2 + 3r 2 h ) 
L 2 (r 2 +4a'k) 



(r - r h ) 2 + a(r) =a h + 



2p 2 rf l a h (pi 



3(r 2 +Aa'k)(2p-l 



-(r - r h 



,2p-l 



+ 



4>{r) = (f>h + 4>i{r - r h f + ... 



(4.3) 



where 



V 



-1 + 



V 



\ 



1 — 2>a\a 2 - 



(k 



ia'k ~» 



Sr 2 



I 



Similar to the nonextremal case, we evaluate the initial conditions (4.3) at r = rh + 10~ 4 
for global tolerance 10 -14 , integrating towards r — > oo. The large r expansion of the extremal 
solutions is still given by the expression (3.17). Similar to the non-extremal case, the value 
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of ah in the horizon data is not relevant in numerics and the only parameter is <f>\. Again, we 
did not notice the existence of local extrema of the functions N(r), cr(r), (j){r). 

In Figure 4 we have ploted the profiles of a typical extremal black hole with non-zero GB 
term together with the corresponding solution with a' = 0. One can see that a non-zero a' 
leads to a deformation of all metric functions at all scales, which holds also for the dilaton (p. 
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Figure 3: The attractor behaviour is shown for k = 1 extremal black holes with A = —1, a' = 0.1, 
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It is also possible to write a simple generalization of the Bertotti- Robinson solution (2.38), 
with the same line element and the matter fields scalar field 

4> = log(-^f), W 1 = e~ ai ' t>0 q 1 r, W 2 = e~ a ^q 2 r, (4.4) 

«2 - ol\ a\qr, 

the size of £§ being fixed by the cosmological constant and the parameters (71,(72 as solution 
of the equation 

^ + A = e-^Ni 2 + e- a2 N 2 2 . (4.5) 

5. Discussion 

In this paper, we have investigated the construction of black hole solutions in higher derivative 
AdS gravity. This is a self contained paper and we hope that our unified treatment of (non- 
extremal AdS black hole solutions with GB term is useful to the reader. 

In the presence of higher derivative terms the area law is modified. Wald proposed a 
new formula for the themodynamical entropy such that the first law of black hole mechanics 
remains valid. In this paper, we have taken a slightly different route in deriving the first law 
and studied the thermodynamical properties of these black holes. The main tool that we have 
used in the non-extremal case is the counterterm method. We have explicitly constructed the 
counterterms that regularize the action and the stress tensor and applied the holographic 
renormalization method in several concrete examples 10 — we found perfect agreement with 
the results obtained by Wald formalism. The method, as developed in [37], constructs unam- 
bigously the counterterms that render the action finite on any solution of the field equations. 
In other words, the method does not proceed with postulating a set of counterterms and then 
checking that these work for a class of solutions. Rather one first obtains the behaviour of 
the on-shell action on the most general solution and then from here one finds the set of coun- 
terterms. It would be interesting to understand better how the counterterms we proposed fit 
in this general framework. 

AdS spacetime is geodesically complete, but the light cones flare out in such a way that 
particles can exit from the space — and also information can come into the space — within a 
finite time. Thus, AdS is not a globally hyperbolic spacetime and so the boundary conditions 
should play an important role. Indeed, within the AdS/CFT duality, various deformations of 
the AdS boundary conditions are interpreted as dual to deformations of the CFT. If such a 
CFT exists the theory lies in the landscape of string theory and the bulk theory is manifestly 
consistent as an effective theory, otherwise the theory is part of the swampland [24] . 

A 'ground state' is defined as a state that extremizes the Hamiltonian over the class of 
vacuum states which all have a given boundary topology. It is well known that by using differ- 
ent foliations of AdS space one can describe boundaries that have different topologies affording 

10 We also have checked this method for other solutions with non-trivial boundary topology [39], but we hope 
to present these results elsewhere. 
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the study of CFT on different backgrounds. The diffeomorphisms in the bulk are equivalent 
with the conformal transformations in the boundary, and different boundary topologies are 
related by singular conformal transformations. Therefore, for AdS gravity, the correct varia- 
tional problem requires keeping fixed a conformal structure rather than a boundary metric. 
The variational problem for two derivative AdS gravity was discussed in detail in [40]. This 
paper also analyzes Wald's method for asymptotically locally AdS spacetimes, shows in gen- 
erality that the results of this method agree with the holographic results, and establishes the 
1st law of thermodynamics. 

Our results fit in this general framework. For black holes with a boundary topology of 
R x S 3 (H 3 ) we found additional Casimir-type contributions to the energy. That is in accord 
with the expectations from quantum field theory in curved space: for the Casimir effect, the 
global structure is reflected nontrivially in the ground state of the quantum field. 

We have also constructed extremal solutions and investigated their properties. Wald 
formalism was extended by Sen to extremal black holes. The advantage of this method is 
that the higher derivatives terms can be incorporated easily, but the method can not be used 
to determine the properties of the solution away from the horizon. However, in section 4, we 
have constructed numerical solutions that interpolate between the horizon and the boundary. 
Thus, we were able to safely apply the entropy function formalism to study their properties. 11 

Unlike the non-extremal case where the near horizon geometry (and the entropy) depends 
on the boundary values of the moduli, in the extremal case, the near horizon geometry 
is universal and is determined by only the charge parameters. We have also constructed 
numerical solutions for which the near horizon geometry is AdS 2 x H 3 (T 3 ), though we do not 
present the details here. We have found that, in all these cases, the scalar fields are attracted 
to fixed values at the horizon. This does not come as a surprise since it is known that the 
'long throat' of AdS2 is at the basis of attractor mechanism. 

A detailed analysis of the attractor mechanism and interpretations within the AdS/CFT 
duality can be find in [15]. 12 After embedding in string theory, the moduli flow becomes 
a holographic renormalization group (RG) flow. The idea that the IR end-point of a QFT 
RG flow does not depend upon UV details becomes in the holographic context the statement 
that the bulk solution in the near horizon limit does not depend upon the details at the 
boundary (asymptotic values of the moduli). That is a holographic interpretation of attractor 
mechanism [15]. 

When the scalars potential is not a constant, a general analysis of the attractor mecha- 
nism is difficult. First of all, if the boundary values of the moduli are fixed to a minimum 
of the potential it is not clear how 'to fly' to IR horizon where the moduli may get different 
values depending of charges (the existence of extremal solutions in this case is problematic). 
However, if the potential has flat directions it may be possible to perturb along these direc- 

11 Entropy function formalism was applied to black holes in AdS space in [15, 41]. 

12 In this paper we are interested in AdS black holes, for more details on attractor mechanism in flat space 
one can consult the nice reviews [42] and references therein. 
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tions. Therefore, a discussion of the attractor mechanism for a non constant scalars potential 
should be made case by case. 

Although the focus of this paper has been on solutions in higher derivative AdS gravity, 
it will be interesting to develop a similar technique for asymptotically flat solutions. In 
particular, it will be interesting to find a boundary stress tensor analogous to the one for two 
derivative gravity [43]. 
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A. Wald formalism for Gauss-Bonet action 



The action in presence of generalised Gauss-Bonnet term is of the form: 



/ 



R VP TPV 



.16ttG \6ttG ^ MW 

Using Wald formalism, the entropy in presence of this term is given by 



(A.l) 



S = / d 3 xVh \l + 2K 5 aR + K b [3{R - h ij + 2K^{R - 2h ij R {j + h ij h kl R ikjl ) 
4G J n L 

(A.2) 

where h is the induced metric on the boundary and K$ = 16ttG. 
Here, we present a detailed proof of this formula. 
General Wald formula for the entropy for any Lagrangean L is 

dL 



-2ir / d 3 xVh 
Jn 



dRabcd 



(abCcd (A.3) 



where 7i is the bifurcate horizon and e„^ is the binormal to the bifurcation surface, normalized 



such that e^ u e^ v = —2. We can take 



V/z, - ivV^ (A.4) 



where ^ and r] are the null vectors normal to the bifurcate Killing horizon, with £.7] = 1. We 
will take 

« - 1 (A ' 5) 
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which is null at the bifurcate horizon. Then r\ can be 

r l= ol ~ 7T ( A - 6 ) 

9tt dt dr 

Now with all these definitions, we can proceed to compute the entropy using Wald formalism. 
We can write the Einstein-Hilbert term using its symmetries as 

R = \(g ac g bd -g ad g bc )Ra bcd (A.7) 

So the leading piece in entropy is 
R 2 part: 

dR 2 

MM**** = ~ AR (A ' 9) 
So, the contribution to the Entropy is,= 

S 1 = -^ J d 3 xVh2K 5 aR (A.10) 



R ij Rij part: 



djR^Rjj) _ or ,ij„l>-ls"j;l'j;c X ,l, 



dRabcd 



-Cabled = ^R %3 g ^i^l^j e ab^cd 



2R bd € a b€d 



(A.ll) 



(A.12) 

= -2R bd (CbVd + ^b), (A.13) 

where we have used the definition of binornal. Now, using the following relation between the 
induced metric and the original metric 

hbd = gbd ~ {CbVd + ZdVb) ( A - 14 ) 
we get that the contribution to the entropy is 

S2 = IgJ d3x ^ K ^( R ~ habR *b) ( A -!5) 
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R^Rijki part: 



d(R tjkl Rjjkl) omkUa!b-!C!d 

dRabcd e ^~ 2R3 W 3***1 



= R abCd (ZaVb-Z b Va)(ZcVd-ZdVc) 

= -2R abcd {g ac - h ac ){g bd - h bd ) 



(A.16) 
(A.17) 



(A.18) 

= -2(12 - 2h bd R bd + h ac h bd R abcd ). (A.19) 

So the contribution to the entropy is 

S3 = IgJ d3x ^ 2K ^( R ~ 2hM Rbd + h ac h bd R abcd ) (A.20) 

Thus we get the net entropy due to the presence of the GB term in the action as S = 
So + Si + S2 + S3 — we used this expression in section 2. 
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